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ABSTRACT 

Context. CoRoT light curves have an unprecedented photometric quality, having simultaneously a high signal-to-noise 
ratio, a long time span and a nearly continuous duty-cycle. 

Aims. We analyse the light-curves of four bright targets observed in the seismology field and study short-lived small 
spots in solar-like stars. 

Methods. We present a simple spot modeling by iterative analysis. Its ability to extract relevant parameters is ensured 
by implementing relaxation steps to avoid convergence to local minima of the sum of the residuals between observations 
and modeling. The use of Monte-Carlo simulations allows us to estimate the performance of the fits. 
Results. Our starspot modeling gives a representation of the spots on these stars in agreement with other well tested 
methods. Within this framework, parameters such as rigid-body rotation and spot lifetimes seem to be precisely de- 
termined. Then, the lifetime/rotation period ratios are in the range 0.5 - 2, and there is clear evidence for differential 
rotation. 
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1. Introduction 

Analyzing stellar activity is important for many top- 
ics. Stellar activity has been studied to increase the ef- 
ficiency of the detection of Earth-like planetary transits 
(jBonomo fc Lanza 2008[) . Tracers of activity can be moni- 
tored with different techniques. White light photometry is 
not the most sophisticated method, and is affected by many 
well known drawbacks (e.g., poor determination of the stel- 
lar inclination, non- uniqueness of the solution). However, 
this method is being revised by new space data, as pro- 
vided by the MOST satellite (jWalker et al. 2003P and now 
by CoRoT. The space mission CoRoT provides us with 2 
years of low noise, long duration, high duty cycle, continu- 
ous time series. The quality of the data, compared to previ- 
ous results, allows us to investigate the activity of solar-like 
stars, and forces us to revise the performance of white light 
photometry. 



The review paper by Collier Cameron (2002) 



summa- 
rizes the method for observing differential rotation, with 
a clear focus on the development of the Doppler imag- 
ing technique. Results, due to the high capability of 
this instrumental technique, were obtained on rapidly 
rotating stars (jColli er Came ron 2007|) . such as the de- 
termination of the magnetic map of a rapidly rotat- 
ing, very-low-mass, fully convective dwarfs through to- 
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mographic imaging from time series of spectropolarimet- 
ric data (jPonati et al. 2006)) . or observations of slowly- 
rotating but active giants (jAuriere et al. 2008)) . Applying 
tomographic imaging techniques requires a study of active 
stars with a magnetic field strong enough to produce de- 
tectable circularly polarized signatures ()Morin et al. 2008)) . 

Space-borne precise photometry opens new windows, as 
already shown by MOST observations. Differential rota- 
tion has been detected by MOST in the light curve of e 
Eridani, a young Sun-like star hosting a planetary compan- 
ion ()Croll et al. 2006p . and in Ceti ([Walker et al. 2007)) . 
Both targets have magnitude variations of about 0.01. With 
CoRoT, it becomes possible to investigate in great de- 
tail activity in solar-like stars with much smaller varia- 
tions. It is for example possible to examine a star such as 
HD 49933, that clearly shows a signal modulated by rota- 
tion IjMosser et al. 20051 jAppourchaux et al. 2008] ), but no 
spectropolarimetric signature (Catala, private communica- 
tion). 

This work addresses the study of the light-curves of 
bright stars with short-lived spots. In most cases, these 
stars exhibit small spots, comparable to the solar spots 
at maximum activity. We have limited the analysis to the 
bright solar-like targets observed in the seismology field. 
For these targets, we have investigated the possibility of 
deriving precise information on the stellar rotation, includ- 
ing possible differential rotation, spot lifetimes and stellar 
inclination. As was done for the fir st analysis of th e MOST 
data, we use the analytic model of Dorren (1987) to model 
the signature of the spots and investigate the case of small 
variations due to many short-lived spots. 
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Observations and light curves are presented in Sect. [51 
The relevant star and spot parameters for the modeling 
are explored in Sect. [S] The performance we may achieve, 
derived from the analysis of synthetic time series and hare- 
and-hounds exercises, is given in Sect. H) Section [5] is de- 
voted to results on the CoRoT time series and to discus- 
sions, and Sect. [6] to conclusions. 

2. Light curves 

We chose to focus this analysis on the main tar- 
gets of the CoRoT seismology field already observed: 
HP 49933 dAppourchaux et al. 20081), HP 175726 



HD49933 



(jMosser et al. 2b09|) HP 181420 (|Barban et al. 2009P 
and HP 181906 (|Garcia et al. 2009^ . The list of the 
stars and their characteristics are presented in Table 
[TJ We worked with the 'Level 2' light curves, ready for 
the scientific analysis, deli vered by the CoRo T pipeline 
after nominal corrections (|Samadi et al. 20W|) . For spot 
modeling, the 32-s initial sampling of the time scries is 
much too dense, so that we rebin the data to one point per 
CoRoT orbit (about 6184 s). With this sampling time, the 
time coverage provided by CoRoT gives a nearly 100% 
duty cycle, except for very short periods when data were 
lost due to telemetric problems between the satellite and 
the ground antennas. This resampling is coherent with our 
spot modeling : it remains tight enough compared to the 
stellar rotation period and has no influence on the spot 
map since the relevance of the fitting process exposed in 
paragraph [3]3] requires the number of spots to be limited. 

Table[T]gives the value of the duty cycles 7732 s and yyorbit, 
respectively calculated with a rapid or a slow sampling. 
The duty cycle ?7orbit useful for the spot analysis depends 
mainly on the number A'intcr of long interruptions (longer 
than 1 CoRoT orbital period) in the data transfer from the 
satellite to the ground. This number is in all cases less than 
or equal to 3. 

The signal-to-noise ratio gains much with such an un- 
dersampling compared to the initial data. The relative stan- 
dard deviation of the light curves due to the photon noise is 
then about 10 ppm for the considered targets, much below 
the intrinsic stellar variations. 

The relative flux variations were then calculated assum- 
ing a reference level $0 close to the observed maximum stel- 
lar flux. In order to remove the non-stellar drifts related to 
the ageing of the CPP, we had to remove a polynomial 
term, close to the linear trend characterized by d In / dt — 
-5.32 10-5 d-i reported by [Auvergne et al. (2009)[ The 
maximum degree of the correction takes into account the 
part of the variation that is correlated for all light curves 
of a given run. 

We chose to analyze 60-day long time series (except for 
HP175726, observed for 27.2 days during a short run). For 
stars observed during the longer than 120-day runs, we 
had then the possibility to split the series and compare 
completely independent subseries. Figure [1] shows the light 
curves. Except for HP 175726, the typical relative varia- 
tions, of about 2 to 3mmag, are similar to the active Sun 
( jAigrain et al. 2004D . 

3. Modeling 

The form of the time series indicates that the light curve is 
modulated by a large number of small short-lived spots. We 
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Fig. 1. Time series of the relative flux variations of 4 
CoRoT bright targets observed in the seismology field. 
Time sampling has been reduced to 1 point per CoRoT or- 
bit (6184s). Only a 60-day sample is shown for HP 181420 
and HP 181906. 
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Table 1. Targets and time series properties. 



star 


type 


mv 


be^ 


jinning 


length 


'732 s 




J/orbit 








(date) 


(CoRoT day) 


(day) 


(%) 




(%) 


HD 49933 


F5V 


5.78 


31 jan 07 


2587.0 


60.7 


93.1 





100 


HD 175726 


GOV 


6.72 


11 apr 07 


2657.2 


27.2 


90.2 





100 


HD 181420 


F2V 


6.57 


11 may 07 


2687.1 


156.6 


91.7 


3 


99.86 


HD 181906 


F8V 


7.6 


11 may 07 


2687.1 


156.6 


91.4 


3 


99.86 



have chosen to expUcitly address this case , with a simple 
model following the formalism developed by Dorren (1987)[ 
modeling the colder stellar patches by individual circular 
spots. We do not introduce either penumbra or facula. As 
in the data observed with MOST and recently analyzed, 
we have to specifically solve for the elements of inclination, 
differential rotation and equatorial rotation period. 

The parameters of interest are, for the star, its inclina- 
tion angle j, the unspotted stellar flux $0, the limb darken- 
ing parameters and the temperature contrast between the 
spots and the stellar photosphere. A spot is characterized 
by the set of an epoch for the maximum spot intensity, a 
longitude, a latitude, and an angular size. Rotation and life- 
time of the spots have an intermediate status: depending on 
the analysis, they were set to the same value for the whole 
star or left free for each spot. The stellar global parameters 
considered as fixed inputs during each spot modeling were 
tested on extensive grids. A large set of values were tested, 
so that their possible variations were precisely analyzed. 



3.1. Limb darkening; spot temperature 

The linear limb darkening function introduced in the model 
is derived from Diaz-Cordoves et al. (1995) according to 
each stellar type. We do not study the influence of the limb 
darkening parameters on the fltting process, since the prob- 
lem in white light is highly degenerate. Due to the fact that 
the spots for the considered set of stars are small, with a 
radius of a few degrees or less, they behave essentially as 
points. Then, a larger spot radius will compensate for an 
underestimated spot temperature deficit, and vice versa. 
Moreover, white light photometry is not suitable for a pre- 
cise determination of those parameters. 



3.2. Unspotted flux, unspotted area 

The value of the unspotted stellar flux is an important issue 
for this study. A too high unspotted level yields for exam- 
ple overestimated spot latitudes, and inversely. Its deter- 
mination therefore has a large incidence on all parameters 
sensitive to the spot latitude; we verified that the incidence 
on the mean rotation and the spot lifetime is limited. 

There is unfortunately no way to obtain an estimate of 
this parameter, and it is not possible to disentangle it from 
the flux variations due to the aging of the CCD camera or 
due to the perturbations caused by high-energy particles 
impacting the CCD. Improving this point will require bet- 
ter knowledge of the CoRoT ultimate photometric precision 
for long-term variation. This work will be undertaken in a 
further work dedicated to differential rotation. 



3.3. Lifetime and number of spots 

First simulations have verified that, as in the solar case, the 
lifetimes of the small stellar spots are short, and comparable 
to the stellar rotation period. We chose to model the time 
evolution of the spot contrast C according to a Gaussian 
law: 



C{t) = Cmax exp 



- In2 



t - tr 



(1) 



with imax the epoch at maximum contrast Cmax, and r the 
lifetime. All spot contributions construct the flux variation: 



(2) 



where $0 is the unspotted flux introduced in Section[51 The 
total number of spots TVs during a lapse of time V has been 
chosen so that there are on average n spots per rotation. 
This is ensured by the relation: 



(3) 



The spot density n is about 2 to 3 per rotation period. A 
larger number of spots reduces the residuals of the best fit, 
but degrades its significance, by introducing, for example, 
large groups of close spots, whose parameters are highly 
degenerate, hence imprecise. 



3.4. Stellar inclination and spot latitudes 

The information on latitude derives from the fact that the 
transit time of one spot depends mainly on this parameter. 
Assuming that, due to possible degeneracies, the time of 
visibility is the only source of uncertainty, an uncertainty 
6i on the inclination i translates directly into an uncertainty 
5\ on the latitude A according to the conservation of the 
ratio: 



tan A tan(A + (5A) 
tan i tan(i -t- 5i) 



(4) 



As a result, underestimating the inclination leads to an un- 
derestimation of the spot latitudes, and vice versa. Even if 
the way the fitting process determines the latitude is more 
complex due to possible interactions between close spots, 
relation (|4]) is useful in estimating the uncertainty: a 10° un- 
certainty on the star inclination may translate into a much 
higher uncertainty on the spot latitudes, which may ham- 
per any precise determination of the differential rotation 
rate of the star. 
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Table 2. Different fitting methods tested, with fixed or free 
parameters. Free parameters are fitted; fixed parameters 
are tested on a grid of values. 



rotation is fixed by 2 parameters in method D, according 
to the rotational law: 



r(A) = 



Fitting method 


A 


B 


C 


D 


inclination i 


fixed 


fixed 


fixed 


fixed 


number of spots TVs 


fixed 


fixed 


fixed 


fixed 


lifetime 


free 


fixed 


fixed 


fixed 


rotation 


free 


free 


fixed 


fixed 


differential rotation 









fixed 


spots 


free 


free 


free 


free 


stellar parameters 


2 


3 


4 


5 


spot parameters 


6 


5 


4 


4 


parameters 




3+5iVs 




4+4iVe 



1 - K sin^ A 



(5) 



4. Fit and performance: hare-and-hounds exercises 

Fitting was performed by an iterative linear search first, 
followed by a Powell search ([Fletcher fc Powell 1963|) . the 

function to be minimized being the sum of the squared 
residuals. Measurements uncertainties being nearly con- 
stant due to the high photometric stability provided by 
CoRoT, the reduced are simply obtained by dividing 
the value by a constant term corresponding to the mea- 
surement uncertainty. As usual in photometric spot mod- 
eling, these reduced are high. During an iteration step, 
all spot parameters have their values fixed by the preceding 
step, but the spot whose parameters are fitted. The extent 
of the exploration of these parameters decreases during the 
iterative process. 

Solutions are relaxed to prevent the solution from 
stalling in a local minimum and to quantify the non- 
uniqueness problem. During the first steps of the iteration, 
relaxation consists of resetting the spot latitudes to zero or 
in adding large perturbations on all spots parameters. In 
the following steps, we control the extent of the exploration 
of the parameters in order to simultaneously enhance the 
precision and reach the convergence as rapidly as possible. 
Globally, changes become increasingly smaller, but simi- 
larly to simulated annealing, the control includes random 
changes. These relaxation steps are achieved automatically 
during the fitting process, according to a global empirical 
scheme optimized when treating the hare-and-hounds exer- 
cises. This scheme defines the number of iterations with a 
given step, the criterion for allowing the next step, the de- 
crease of the changes of the exploration, and the occurrence 
of the relaxation events. During the fitting process, spots 
can be removed when their contribution becomes negligi- 
ble, due to a too small angular radius and/or a too high 
latitude. 

With 4, 5 or 6 free parameters for each spot, depend- 
ing on whether the rotation and lifetime are fixed or not 
(Table ^ , and typically 50 spots for each 60-day long se- 
quence, the number of the parameters to be fitted varies 
from about 200 to 300. Fitting method A, with both spot 
lifetimes and spot rotations considered as free parameters, 
was used as a first step. This method A is useful for de- 
termining the relevant ranges of these parameters, hence 
for defining the grid of the values tested by the other fit- 
ting methods. It shows that spot lifetimes are mainly dis- 
tributed in a restricted range. Therefore, the lifetime t was 
then considered as a fixed stellar parameter for all other fit- 
ting processes. Spot rotation period is free in fitting method 
B, and fixed to a single value T in method C. Differential 



with Teq the equatorial rotation period and K the differen- 
tial rotation rate. 

Before applying the spot modeling to the CoRoT tar- 
gets, we have tested it on synthetic light curves. An exten- 
sive study has been provided, with blind analysis. Synthetic 
light curves were constructed by one of us and given to the 
analyzer without any information. The parameters of the 
synthetic light curves were varied in order to mimic the va- 
riety of the parameters estimated in the real stars. In order 
to cover a set of cases at least as large as the observations, 
we tested inclinations in the range [10 - 90°], lifetimes r 
with ratios r/T in the range [0.5-2], and mean spot radius 
below 0.5°. 

The method proved to be robust for retrieving the stel- 
lar parameters according to the results of these hare-and- 
hounds exercises. With numerous tests, we were able to de- 
rive realistic uncertainties, derived from the standard devi- 
ation of the difference between the input and output values. 
These uncertainties on the rotation, spot lifetime and star 
inclination are presented in the following paragraphs. 

4.1. Rigid-body or differential rotation 

Simulations have shown that the determination of the mean 
rotation T is robust. The uncertainty we get on this pa- 
rameter is about 0.5 % in most cases, up to 1 % in the un- 
favourable case of very short-lived spots or of stellar incli- 
nation close to the pole-on case. Of course,_when differential 
rotation is present, the determination of T, defined as the 
mean value of the spot rotations, is possibly affected by the 
distribution in latitude of the spots. Simulations show that 
T remains however a relevant parameter. 

The determination of the differential rotation parameter 
K is not straightforward. When the star is seen pole-on or 
edge-on, the degeneracy on the latitudes precludes any reli- 
able determination of the coefficient K. Values of about 0, 
15 an 30 % were tested. Better determinations are found for 
mid-inclination models. The determination of K is usually 
underestimated. Precision in the hare-and-hounds exercises 
is limited to about 70%. A large part of the error comes 
from the uncertainty on the star inclination. 

4.2. Spot lifetime 

As for the rotation, we estimate the performance of the 
method when trying to retrieve the spot lifetime. The life- 
time appears to be slightly overestimated if too few spots 
were included in the modeling, and underestimated in the 
opposite case. When the spot lifetime is long compared to 
rotation, the fitting process is less precise, due to inter- 
actions between spots. These interactions may introduce 
uncontrollable degrees of freedom. 

The estimate of the mean spot lifetime r in the syn- 
thetic time series is robust, and not significantly affected 
by uncertainties on other parameters. The precision we get 
on this parameter is about 10 %, in the case where the ratio 
r/r is in the range [0.5 - 1.5], and about 20% for larger ra- 
tios. When lifetimes become much shorter than the rotation 
period, the error bars on all parameters increase. 
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Table 3. Output parameters derived from the spot modeling. 
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star 



inclination 



spot 



(deg) 



Spot contrast 

(Cm ax) 

(7oo) 



spot angular radius 
median range 
(deg) 



lifetime 



mean rotation 

r 

(day) 



HD 49933 
HD 175726 
HD 181420 
HD 181906 



50 ± 25° 
55 ± 25° 
60 ± 25° 
45 ± 25° 



55 ±10° (a) 



45 ±4' 
24 ±3' 



(b) 
(c) 



1 

8 

0.8 
1 



1.8 
4.7 
1.6 
1.8 



1 ^ 
3 - 
0.8 
1.0 - 



2.5° 

» 7° 
->2° 
> 2.5° 



2.5 - 
5 - 
1.5 
1.5 



3.5 

7 

^ 3 
^ 3 



3.45: 
3.95^ 



-0.05 
^+0.1 
-0.1 
9 9t-+0.03 
^■^'-'-0.01 
9 7-1 +0.03 
^- ' -'--0.01 



Inclinations derived from the spot model are compared to the seismic values inferred from the fit of the Fourier spectra of p modes: 
(a) Appourchaux et al. (2008); (b) Barban et al. (2009); (c) Garcia et al. (2009). 
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Fig. 5. Low-frequency spectrum of HD181906. The rota- 
tional signature exhibits 2 peaks at about 4.0 and 4.4 /iHz 
(2.9 and 2.6 d, respectively). 



T (day) 

Fig. 6. Histogram of the distribution of the spot rotations 
in HD181906. 



4.3. Star inclination and spot latitudes 

Determining inclinations is well known to be highly dif- 
ficult with white light photometry. The hare-and-hounds 
exercises show that it is however possible to recover cor- 
rected information on the latitude for intermediate lati- 
tudes. Determining the latitude for pole-on on edge-on cases 
is however very imprecise, due to degeneracy. 

To fit the inclination, different values were exhaustively 
tested between 20 and 90°. The precision on the latitude 
in the synthetic time series is about ±10° for inclinations 
in the range [30 - 60°], increasing to ±20° in the range [70- 
90°]. When the star is seen edge-on, no information can be 
extracted to derive the spot latitude, since all spots transit 
the star in half a rotation period. When the star is seen 
pole-on, there is no further transit, and smooth modula- 
tions are difhcult to fit. It is difficult to distinguish a quiet 
star from an active one seen pole-on. In all cases, the de- 
termined stellar inclination angle depends significantly on 
the assumptions of the model. 



5. Results and discussion 

Hare-and-hound exercises presented above make us confi- 
dent to analyze the CoRoT time series. Spot modeling has 
been applied to 60-day long times series of the CoRoT as- 
teroseismic targets (Fig. [21 [21 S]) . The results we obtain for 
the star inclination, the spot size distribution, the lifetime 
of the spots and the mean rotation period are presented 
in Table [31 and discussed in this Section. Deriving reliable 
error bars being not straightforward in the case of a large 
number of parameters, we applied three different methods. 



- The tests made on synthetic time series allow us to quan- 
tify the goodness of the fit and to extract confidence inter- 
vals from the chi-square fitting. Assuming Gaussian statis- 
tics, due to the large amount of data, the relative likelihood 



of a set S of parameters can be expressed as Croll (2006) 



ps cx exp 



(xl 



Amin 



2T 



(6) 



with Xmin the minimum value reached for the best fit, and T 
the virtual temperature as used in the case of Markov chain 
Monte Carlo methods. This virtual temperature should be 
1 when the value of the reduced chi-square is about 1. 
With large values of Xmin' ^^'^ range fro m 50 to 100 
in most cases, we should use T — x'^^^\n in [Croll (2006) 



However, hare-and-hounds exercises show that the 1-ct un- 
certainty derived from Eq. [51 (at 68 % of the maximum value 
reached by ps) is then significantly exaggerated. In fact, an 
irreducible contribution to the values comes from high- 
frequency signatures that the modeling cannot reproduce; 
it would require the introduction of too many small spots, 
that would translate into quantitatively good fits, but qual- 
itatively poor, or the introduction of hot faculae. In order 
to reconcile the performance of the method derived from 
the hare-and-hounds exercises and Eq. [51 we have empiri- 
cally considered T ~ 0.4Xmin- We are aware that this sup- 
poses that the models tested during the exercises were close 
enough to reality to give reliable indications of the preci- 
sion. 

- A second estimate of the confidence interval was achieved 
by comparing results obtained from tens of fits for each in- 
put configuration, derived from different initial conditions 
and different relaxation tools. 
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- A third estimate was derived from the comparison of in- 
dependent time subseries, available for HD 181420 and HD 
181906. We have verified that the independent best fit val- 
ues we obtain agree. 

Observed error bars were defined in order to encompass 
all results derived from the fits based on methods B and C, 
method A being used to define the grid of the parameters 
to be tested. These error bars are slightly larger than the 
ones derived from synthetic time series. The partial results 
of method D are consistent with other methods, but need to 
be completed by more fits before being reliably presented. 
This will be done in a forthcoming work, and compared to 
the results given by the maximum entropy method (Lanza 
et al. 2009a). 

5.1. Rotation 

The mean star rotation T can be precisely inferred from the 
fits. The relative precision is of about 1 % (but 2.5 % for HD 
175726, observed for 27.2 d), comparable to the precision 
that can be achieved from the study of the low-frequency 
spectrum (|Baudin et al. 20091) . The presence of differential 
rotation is strong ly suspected in F dwarfs, as reported by 
Reiners & Schmitt (2003), and, in fact, three out of these 
four targets show indications of non-solid rotation. 

HD 181906 represents the clearest example: as already 
noticed (jGarcia et al. 2009|, its low- frequency spectrum 
(Fig. [SI) exhibit 2 peaks, at about 4.0 and 4.4 /iHz. The 
spot modeling also shows 2 peaks in the histogram of the 
spot rotations derived from the fitting method B. The 2 
peaks at respectively 2.6 and 2.9 d correspond clearly to 
the signatures in the Fourier spectrum (Fig. [S]) . They may 
be due to a combination of the rotation period and short 
spot lifetime. We checked that they are more likely due to 
differential rotation, since the latitudes of the spots are cor- 
related to the rotation rate: the mean latitudes of the spots, 
sorted according to the rotation period, increase with in- 
creasing rotation period (Fig. [7]). HD 175726 presents sim- 
ilar features (Fig.[8|). 

We can derive the order of magnitude of the differential 
rotation rate, but with large error bars: K is about 0.25 
for HD181906, and about 0.40 for HD1 75726, in agre ement 
with expected values for such targets (jReiners 2006p . The 
case of HD 49933 will be exhaustively considered in a forth- 
coming work. 



5.2. Spot lifetime 

If our spot model is correct, then as expected from the pre- 
liminary synthetic analysis, the determination of the spot 
lifetimes gives convergent, hence reliable results. The values 
we obtain fit in the range 0.5— s-2 times the rotation period. 
The lifetimes of the spot may increase with increasing ro- 
tation period (Fig. [9]) and with increasing spot size. 

The uncertainties we derive on the spot lifetimes vary 
from ±20 to ±50%. Different reasons may explain these 
large values: the ad-hoc model to describe the evolution 
with time (Eq. [1]), and dependence of the hfetime on the 
spot si ze. The study of the low-frequency spectrum of HD 
49933 C Baudin et al. 2009!) gives similar results, which en- 
hances the confidence we may have in the description of the 
spot map provided by the m odel. We also applied the max- 
imum entropy spot model of Lanza et al. (2009a) to derive 
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llotltudel (deg) 

Fig. 7. Rotation-latitude relation, for HD181906. The size 
of the solid black or grey squares (for spots respectively 
located in the northern or southern hemisphere) is propor- 
tional to the spot size. The solid line indicates the best fit 
of differential rotation according to Eq. O 
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Fig. 8. Same as Fig. H but for HD175726. 



information on the active region evolution of HD49933. The 
time series were subdivided into subsets of one rotation pe- 
riod (i.e., 3.4 days), each of which was fitted with a fixed 
distribution of spots, optimized to fulfil the maximum en- 
tropy regularization criterion. Therefore, this approach is 
sensitive only to evolutionary timescales longer than one 
rotation period, in contrast to the discrete spot modeling 
approach, which can resolve shorter spot lifetime. In spite 
of such a difference, the results obtained with the two meth- 
ods are comparable. In the ME model of HD49933, most 
of the spots have lifetimes between 5 and 7 days, with only 
10 — 15 % showing a lifetime up to 7 — 10 days during the 
short run time series. 



5.3. Inclination 

Precisely determining the star inclination has proven to be 
difficult, and we cannot exclude that the results are model 
dependent. We did not use any indication derived from the 
stellar v sinz. Error bars for all targets are at least of ±25°. 
This means that the fitting process is much less efficient 
than for the rotation period and spot lifetime. Surprisingly, 
all four values we find (Table [3]) correspond to inclinations 
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Table 4. Stellar inclinations. 

star wsini R/Rq 2nR sin ispot /v sin i T !{ V sin 2, T) ispot 



(km s-^) (d) (°) 



HD 49933 


9.5±0.3 (a) 


1.39±0.02 (d) 


5.7 


3.451-- 


27±2 


50±25 


HD 175726 


13.5±0.5 (b) 


l.OiO.l (c) 


3.2 


3.95ir, 


^"-30 


55±25 


HD 181420 


18±1 (c) 


1.6±0.1 (c) 


3.9 


9 9r+0.03 
^•^^-0.01 


30±4 


60±25 


HD 181906 


10±1 (c) 


1.4±0.1 (c) 


5.0 


9 71 +0.03 
^- ' J^-0.01 


23±4 


45±25 



References: (a) Mosser et al. (2005); (b) Kovtyukh et al. (2004); (c) Bruntt et al. (2009); (d) Goupil et al. (2009) 



2.0 2.5 3.0 3.5 4.0 

Mean r'otcit.ion (day) 

Fig. 9. Spot lifetime as a function of the mean rotation 
period. 

in the range [45 - 60°]. According to the preliminary work 
on synthetic data, we suspect a bias in favor of such mid- 
values, but the inclination values derived from asteroseis- 
mology then present the same behavior, or the same bias. 
Such a bias could be due to the fact that rotational multi- 
plets in the case of mid-inclination are more complex com- 
pared to the degenerate cases i = 0° (only £ = modes 
are observable) or 90° (only even values of \£ + m\ present 
non-zero visibility). However, the sample of stars presently 
is too limited to draw conclusions. 

Obtaining reliable measurements on differential rotation 
should be supported by other measurements, by avoiding 
uncertainties due to an inaccurate determination of the stel- 
lar inclination. The best unbiased estimate of stellar incli- 
nation should come from the combined determinations of 
the rotation period T from spot modeling, of the stellar ra- 
dius R determined from the seismic modeling, and of v sin i 
derived from ground-based high-resolution spectroscopy. In 
Table 01 we present the corrected inclinations we obtain ac- 
cording to this method. Up to now, only HD49933 bene- 
fits from precise radius measurement derived from seismic 
modeling. The discrepancies in the inclination are high and 
deserve further work. We should take these new values of 
the inclination to reestimate the differential rotation rate. 
However, as mentioned above, further work is needed to 
obtain reliable estimates and uncertainties. 



6. Conclusion 

The new data provided by CoRoT renew spot modeling 
with white light photometric data. With long duration, high 
duty-cycle, very high signal-to-noise ratio data, we can now 
address the case of short-lived spots in solar-like stars with 
solar-like activity. The only differences between the active 



Sun and this sample of stars are for the types (F2-^G0 
compared to the solar type G2) and the rotation periods 
(2.7^3.9 d compared to a mean value of 27.3d). 

The work presented here is a first step in addressing 
the study of such stars. The metho d we use is sim pler than 
many previous studies, for example Croll (2006) analyzing 
MOST observations, but here, we deal with many short- 
lived modes, which yields a much larger number of free 
parameters. 

The maximum entropy method proposed by 



Lanza et al. (2009a) is capable of fitting the light curve 
with high accuracy because it uses a continuous distri- 
bution of spots, but it cannot follow the evolution of the 
active regions on timescales shorter than a single rotation 
period. Moreover, stellar rotation period and inclination 
of the rotation axis must be fixed in advance to ensure 
that the regularization process converges toward a unique 
and stable spot configuration. In view of our limited 
knowledge of these stellar parameters, we decided to use 
the approach described in Sect. 3 to obtain information on 
stellar rotation rate and inclination, as well as on the spot 
lifetimes, because they arc shorter than or comparable to 
the rotation period in most of the cases. The application of 
the maximum entropy method is foreseen for future studies 
to investigate the evolution of the complexes of activities 
and derive a lower limit for the surface differential rotation 
in our targets, as shown by, e.g., Lanza et al. (2009a, 
2009b). 

The main result for the four CoRoT targets we stud- 
ied is our reliable and precise measurement of the spot 
lifetimes and mean rotation. Of course, this requires that 
our starspot modeling gives an accurate representation 
of the spots on these stars, which seems to be con- 
firmed by the agreement with the maximum entropy re- 
sults and with the analysis of the low frequency spec- 
trum (I Baudin et al. 2009]) . We also find evidence of dif- 
ferential rotation, and show that independent measure- 
ments of the star inclination will be of great help in in- 
ferring the rate of differential rotation. Spot modeling 
may then achieve complementary results to asteroseis- 
mology (Gizon fc Solanki 2004[ |Ballot et al. 2006p . A dedi- 
cated analysis of differential rotation will be presented in a 
forthcoming work. We will also explore other ways to per- 
form the fit of light curves with short-lived spots, such as 
Markov Chain Monte Carlo methods. 

Acknowledgements. We thank the referees, Gordon Walker and Bryce 
Croll, for their advised comments that help improving the paper. 
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Fig. 2. Spot modeling of HD49933 and HD 175726. The dots in each upper panel represent the data binned every CoRoT 
orbit (6184 s), and the solid curve is the best fit model. Dots in each bottom panel represent residuals. Vertical grey lines 
indicate the mean rotation period. 
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Fig. 3. Same as Fig.d but for HD181420 between CoRoT days 2690 and 2750, and between CoRoT days 2770 and 2830. 
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Fig. 4. Spot modeling of HD181906, same as Fig. [31 



